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ABSTRACT 

The  metal  vapor  from  a  wire  exploded  in  a  gas  is  assumed  to  be 
contained  between  the  walls  of  a  hollow  cylinder  which  has  boundaries 
that  are  linear  functions  of  time.  It  is  shown  that  the  presence  of 
the  ambient  gas  can  be  ignored  in  the  qualitative  descriptions  of  the 
flow  determined  from  interf erograms .  Studies  of  the  single  and 
multiple  fringe  interferograms  expected  from  such  a  cylindrical  flow 
are  made.  Graphical  plotting  techniques  are  developed  to  obtain  these 
expected  interferograms,  and  supplementary  mathematical  relations 
describing  the  interferograms  are  demonstrated.  Finally,  comparison 
of  interferograms  obtained  by  these  methods  is  made  with  those  inter¬ 
ferograms  obtained  during  actual  exploding  wire  experiments. 
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I.  INTRODUCTION 


Single  and  multiple  fringe  interferograms  have  been  in  use  in 
recent  years  to  obtain  gas  and  electron  distributions  about  an 
exploding  wire-1"*.  Until  recently  not  much  quantitative  evidence  has 
been  obtained  about  the  flow  of  the  metal  in  these  exploding  wires. 

The  purpose  of  this  paper  is  to  establish  some  geometrical  techniques 
and  mathematical  relations  upon  which  subsequent  analyses  of  this  metax 
flow,  seen  in  the  existing  interferograms,  may  be  based.  This  paper 
offers  certain  flow  hypotheses,  shows  the  resulting  fringe  contours  for 
both  the  single  and  multiple  fringe  cases,  and  compares  these  fringe 
contours  with  some  actual  interferograms. 

II.  FLOW  HYPOTHESES 

A.  Outer  Zone  Cases g 

Figure  1  represents  the  expansion  as  a  function  of  time  of  the 
boundaries  of  the  metal  vapor  from  a  wire  exploding  in  the  presence  of 
a  low  pressure  gas.  We  assume  this  flow  to  be  cylindrically  symmetric 
about  the  wire  axis.  Gas  density,  p  ,  is  assumed  constant  for  r  >  r^. 
In  the  shocked  ambient  gas  such  that  r™  ^  r  >  r^T  the  approximate, 

constant  value  is  assumed.  For  r^  ^  r  ^  rp  the  metal  vapor  is 

represented  by  constant  density  p^  ;  and  for  the  innermost  region, 

where  rp  >  r  ^  0  ,  density  is  taken  to  be  zero.  Note  that  all  of 

the  boundaries  are  assumed  to  be  linear  and  intersect  the  r  axis  at 
time  zero.  This  approximation  fits  the  observed  flow  boundaries  fairly 
closely  and  presents  no  serious  problems  except  near  t  =  0.  Thus  pp 
can  be  represented  approximately  as  a  hollow  cylinder  of  metal  vapor 
of  inner  radius  rp  =  az  and  outer  radius  r^  =  bz.  Axial  distance  is 
related  to  time  through  the  equation  z  ^  2  £co t,  to  being  the  angular 
velocity  of  a  mirror  making  the  streak  photograph  of  the  event  repre¬ 
sented  in  Figure  1,  I  is  the  optical  lever  arm,  and  t  is  the  time 

2 

elapsed  since  the  initiation  of  the  explosion  .  If  mass  is  conserved, 
^References  are  listed  on  page  30  . 


9 


and  the  original  radius  and  density  of  the  wire  are  represented  by  r 

w 

and  p  respectively,  then  pTI  =  p  r.2/(b2-a2)  z2  ,  for  r  4  r.  The  ex- 

W  "  it  W  w  '  '  W 

press ion  for  a  fringe  shift  6  at  radius  r^ ,  for  rays  which  are  normal 
to  the  r  -  z  plane  of  a  cylindrically  symmetric  disturbance  of  outer 
radius  r0 ,  is  given  by  3 


o  1°  (n-n  )sds 

h\  - ^ 


r.  (s2-r.2) 

i  i 


(1) 


where  s  is  the  radial  distance  to  a  point  on  the  ray  passing  through 
the  disturbance  at  r.,  n  is  the  index  of  refraction  at  s,  nQ  is  the 
index  of  refraction  of  the  medium  in  the  reference  beam3,  and  X*  is  the 
vacuum  wavelength  of  the  light  used.  For  the  disturbance  assumed  in 
Figure  1,  n  is  not  a  function  of  s,  and  the  integration  may  easily  be 
carried  out  between  the  limits  indicated  in  Figure  1,  with  the  result 


If  6<ri’z)  =  Si  )Z5  -  P0)[(^3-^)2  -  (a2 -a2 ) 


.1-1 

2 


+  K  p  [3(c2  -o'2  )2  -  3 (b2  -a2  ) 2  -  (a2  -c*2  )2  ] 

o  O 


(2) 


where  in  each  zone  the  Dale -Gladstone  approximation  relating  density 
and  index  of  refraction  has  been  made.  The  Dale-Galdstone  constant 
used  for  the  vacuum  and  ambient  compressed  gas  zones  is  and  is 
used  for  the  metal  vapor  zone.  The  parameter  a  =  r/z  has  been  intro¬ 
duced  for  convenience.  Assuming  that  pQ  is  the  density  of  l/lo  atm 
argon ,  that  Kg  =  K^,  and  that  the  values  for  a,  b,  and  c  are  as  shown 
in  Figure  2,  one  may  plot  the  (r,z)  contour  for  a  given  fringe  shift, 
e.g.,  6=6,  by  solving  Eq.  (2)  for  z,  calculating  z  as  a  function  of 
o^s  of  interest  (c  >  a  ^  0),  and  recalling  r  =  cvz.  This  procedure 
results  in  the  solid  curves  of  Figure  2  for  an  r^  of  5  mils.  The  inner 
solid  curve  represents  the  negative  root  of  z,  the  outer  the  positive 
root. 
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B .  Approximation  to  Outer  Zone  Case. 

Because  p  r2/(b2-a2)z2  p  and  the  bracketed  radicals  are  all  of 
w  '  o 

the  same  order  of  magnitude  for  the  experiments  performed,  we  might 
reasonably  assume  that  the  effect  on  6  of  the  outer  zone  and  inner 
vacuum  zone  is  small.  With  this  assumption,  Eq.  (2)  may  be  written 


6 


2KMrwX 
X*(b2  -a2)z 


[(b2  -o2  )2 


(a2  -o'2  )  ]2  . 


(3) 


Following  the  same  plotting  procedure  and  using  the  same  numerical 
fringe  shift  number  as  before,  Eq.  (3)  yields  the  dashed  curve  in 
Figure  2.  The  adjacent  solid  curve  is  quite  similar  in  shape  to  the 
dashed  curve,  differing  only  in  magnitude.  Defining  the  difference 
between  ordinate  values  for  the  two  lower  curves  as  Ar,  it  can  be  shown 
that  at  a  fixed  z  as  6  ->  «,  Ar  mo  no  tonic  ally  approaches  0*  Because  the 
6  chosen  in  the  figure  (6  =  6)  is  small  compared  to  those  likely  to 
occur  in  actual  experiments,  the  difference  between  the  lower  curves  is 
about  the  largest  that  might  be  expected.  Since  this  difference  is 
reasonably  small,  the  use  of  Eq.  (3)  in  place  of  Eq.  (2)  is  assumed  to 
be  justified  in  representing  the  lower  curves.  The  (r,z)  contour  from 
Eq.  (3)  does  not  have  an  outer  curve  to  compare  with  that  from  Eq.  (2). 
In  any  case  the  outer  curve  does  not  occur  at  z  values  of  current 
interest  (z  <  5  cm),  nor  has  it  been  observed  experimentally.  Conse¬ 
quently,  the  fact  that  Eq.  (3)  yields  no  outer  curve  is  not  considered 
a  serious  shortcoming  in  its  use  in  preference  to  Eq.  (2).  We  will  use 
Eq.  (3)  therefore  in  subsequent  discussions. 


III.  GEOMETRICAL  CONSTRUCTION 
A.  The  Function  G(o') 

Defining  the  bracketed  term  of  Eq.  (3)  as  G(cr)>  we  recall  how  the 
function  is  generated.  In  Figure  3  for  0  £  a  £  a,  G(a)  is  represented 
as  the  length  of  the  line  cd,  which  is  the  difference  in  ordinate 
values  for  two  circles  of  radius  a  and  b  and  abscissa  a.  For  a  <  cy  ^  b, 
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G(a)  can  be  represented  as  the  ordinate  distance  to  the  radius  b  circle. 
Hence 


G(a)  =  (b2-*2)2  -  (a2-o*)2, 


=  (b2-^)2. 


0  s  a  ^  a 

(4) 

a  ^  or  £  b . 


G(a)  versus  a  is  also  shown  in  Figure  3  as  the  heavy  curve.  Note  that 

G  =  (b2-a2)2,  G(o)  =  b  -  a,  and  G(b)  =  0.  One  can  easily  verify  that 
max 

G'(a)  =  +  oo,  G'(b)  =  -  oo?  and  G'(0)  =  0  as  shown  in  the  figure. 

B.  Single  Fringe  Construction. 

From  Eqs.  (4),  Eq.  (3)  may  be  written  as 


z  =  KG(c*)/6 


(5) 


where  the  definition  of  K  is  evident  by  inspection.  Keeping  6  fixed, 
one  may  plot  the  lines  r  =  az  and  z  =  KG(o')/6  for  all  a 1  s  of  interest. 
Where  these  curves  intersect  for  the  same  o'  is  a  point  in  the  r  -  z 
plane  for  the  chosen  fringe  shift  6.  The  locus  of  intersection  points 
therefore  describes  the  fringe  shift  contour  in  the  r  -  z  plane.  Such 
a  geometrical  procedure  is  followed  in  Figure  4,  where  the  bottom 
abscissa  represents  the  quantity  KG(cv)/6.  The  fringe  shift  contour  is 
traced  out  in  the  top  half  of  ' the  figure.  Fringe  shift  patterns  of  this 
type  are  frequently  observed  experimentally.  Figure  5  is  a  single 
fringe  interferogram  of  a  copper  wire  exploded  in  l/8  atm.  argon.  The 
wing-like  structure  of  Figure  4  is  apparent  in  the  region  of  Figure  5 
close  to  the  r  =  0  axis. 

Interfere grams  of  wires  exploded  in  a  vacuum  indicate  the  absence 
of  the  inner  boundary  defined  by  the  line  r  =  az,  i.e.,  the  expansion 
takes  the  form  of  a  solid  cylinder  with  a  radius  that  is  increasing 
linearly  with  time.  Figure  6  is  an  interferogram  of  this  vacuum 
explosion.  It  should  be  compared  with  Figure  7  which  is  a  plot  of 
Eq.  (3)  for  a  =  0,  b  =  1. 
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C .  Multiple  Fringe  Interf erograms . 


The  fringe  shift  6  is  given  by  6  =  r/W  -  rQ/W  where  r  is  the 
radius  of  the  fringe  in  the  disturbed  region,  and  rQ  is  the  radius  of 
the  same  fringe  in  the  undisturbed  region.  Fringes  in  the  undisturbed 
region  are  aligned  parallel  to  the  z  axis  and  equally  spaced  a  distance 
W  apart.  Then  from  Eq.  (5) 


(r  -  r^)z  =  KWGCa) . 

Dividing  Eq.  (6)  by  |rQ|  and  defining  r  =  r/|r  | 

(r  T  1)2  =  KWG(a-)/r02  . 

The  upper  sign  represents  the  case  where  rQ  >  0, 
Also 


(6) 

and  z  =  z/|rQ|,  yields 

(7) 

the  lower  where  rQ  <  0. 


r  =  ±  xz  , 


(8) 


where  x  =  |c*|  . 

Eq.  (7)  represents  equilateral  hyperbolae  with  origins  at  r  =  ±  1, 
z  =  0.  Note  that  with  such  hyperbolae  the  radial  distance  between  the 
origin  and  the  hyperbola  measured  along  the  (r  T  l)  =  z  line  is  given 
by  (2KWG/r02)^  .  Thus  for  rQ  >  0,  one  may  plot  Eq.  (7)  for  a  given  a, 
r  ,  K  and  W  and  obtain  a  hyperbola  in  the  r  -  z  plane.  One  may  also 
plot  Eq.  (8)  and  the  intersection  point  of  this  straight  line  with  the 
hyperbola  corresponds  to  the  value  of  (r,  z)  for  the  a  chosen.  Select¬ 
ing  all  or's  of  interest,  one  may  thus  determine  the  contour  of  a  fringe 
as  it  passes  through  the  region  of  disturbance.  This  procedure  is 
followed  in  Figure  8.  On  the  inclined  axis  the  radial  distance  is 
plotted  versus  ot.  Figure  9  is  a  multiple  fringe  interferogram  taken  in 
the  presence  of  argon  at  l/l6  atmosphere.  The  hook-like  behavior  of 
the  fringe  in  Figure  8  can  be  seen  in  the  fringes  of  the  upper  portion 
of  Figure  (9) • 
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A  similar  plotting  treatment  for  rQ  <  0  is  carried  out  in  Figure 
10.  In  this  case  there  are  two  branches  of  the  fringe  contour  corre¬ 
sponding  to  two  points  of  intersection  of  Eqs.  (7)  and  (8).  For  some 
values  of  ot  these  two  branches  may  intersect  as  shown  in  Figure  11  to 
form  an  "excluded  region",  i.e.,  a  region  where  for  a  particular  range 
of  a,  no  fringes  can  occur.  Geometrically,  this  results  from  the 
abscissa  values  of  a  set  of  hyperbolae  being  everywhere  greater  than 
those  for  the  corresponding  set  of  r  =  -  xz  lines,  i.e.,  no  inter¬ 
section  occurs.  Evidence  of  this  excluded  region  has  not  been  found 
experimentally . 


IV.  MATHEMATICAL  RELATIONSHIPS 
A.  Single  Fringe  Case. 

From  Eqs.  (4)  and  (5) 


(2c^-bs) 

O' 


0  ■s  a  <  a. 


a  <  or  <  b. 


(9) 


Hence  for  single  fringe  interferograms 


dr 

—  =  a,  b,  +  oo,  at  a  =  a,  b,  0,  respectively. 

CL  2# 

These  slope  values  correspond  qualitatively  to  those  of  the  single 
fringe  interferogram  of  Figure  5- 

B.  Multiple  Fringe  Case. 

From  Eqs.  (7)  and  (8)  one  may  solve  for  z  and  obtain 


(-1)1  +  (-l)j  P 

(-l)k  2x 


lU 


z 


(10) 


where  P  =  1  +  (-1)^  4GWKx/r^2  .  The  values  for  the  possible  fringes, 
seen  in  composite  in  Figure  12,  are  given  in  Table  I. 


Table  I.  Matrix  of  Indices 


Region 

i 

j 

_ _ , 

k 

Fringe 

r  >  0,  ?  >  0 
o  ’ 

mmm 

i 

r  >  0,  ?  >  0 
o  ’ 

even 

odd 

even 

not  allowed 

rQ  >  0,  r  <  0 

even 

even 

odd 

not  allowed 

r  >  0,  r  <  0 
o 

even 

odd 

odd 

not  allowed 

r  <  0,  r  >  0 

0 

odd 

even 

even 

ii 

r  <  0,  ?  >  0 
o 

odd 

odd 

even 

not  allowed 

r  <  0,  r  <  0 
o  * 

odd 

even 

odd 

m 

r  <  n  r  <  0 
o 

odd 

odd 

odd 

IV 

When  k  is  odd,  Eq.  (10)  has  imaginary  roots  when  1  <  4  GWKx/r_2.  This 
means  no  representation  in  the  real  r  -  z  plane.  These  imaginary  roots 
correspond  to  the  "excluded  regions"  of  Figure  11.  Note  that  at  the 
boundaries  of  the  excluded  regions. 

1  =  4GWKx/rQ2  (11) 

so  that  the  upper  and  lower  branch  always  intersect  at  r  =  -  l/2. 

The  slope  dr/dz  can  be  determined  as  in  Eq.  (9).  The  result  is 


where 


dr 

dz 


_ 2x3M _ 

(-l)k2xM  -  P  +  (-l)1+J+1p2 


(12) 


M(x)  =  (KW/ro2)  d(xG)/dx  . 
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The  values  for  i,  j ,  k  for  a  given  fringe  are  again  found  in  Table  I. 
Hence  at  the  intersection  of  the  upper  and  lower  branches  of  the  curve 
below  r  =  0  line,  dr/dz  =  -  xx ,  -  Xg ,  where  xx  and  Xj  are  roots  of 
Eq.  (11 ) .  Also  dr/dz  =  (-l)k|a|,  (-l)k|b|,  at  x  =  |al,  |b|,  respectively 
Also  for  both  curves  II  and  III 


dr 

dzx=o 


~ro  1 
KW  (b-a)  * 


(13) 


Comparing  this  value  with  that  of  the  slope  for  single  fringe 
interferograms ,  we  see  that  the  two  slopes  of  the  single  and  multiple 
fringe  interferograms  can  never  be  equal  at  x  =  0. 

C .  Asymptotic  Behavior. 

The  shape  of  the  fringe  for  large  z  for  the  multiple  fringe  case 
can  be  seen  from  Eqs.  (8)  and  (10).  Combining  these 

fM  -  (it) 

where  i,  j ,  and  k  are  even  or  odd  integers  as  given  in  Table  I,  r(x) 
at  z  =  »  is  given  by  r(o).  For  the  curve  of  Figure  8,  curve  I, 

(r  >  0)  i  and  j  are  even  integers  hence  r(o)  =  1.  That  is,  the  fringe 
after  passing  through  the  region  of  disturbance  will  asymptotically 
approach  its  initial  configuration,  i.e.,  a  line  parallel  to  the  r  =  0 
axis  at  a  distance  of  unity  away  from  the  axis. 

Likewise  for  the  lowermost  curves  of  Figure  10,  curve  TV,  (rQ  <  0) 
i  and  j  are  odd.  Hence  from  Eq.  (l4),  r(o)  =  -  1.  The  lowermost  curve 
thus  approaches  its  initial  position  with  increasing  z.  This  asymp¬ 
totic  behavior  for  rQ  ^  0  is  discernible  in  the  multiple  fringe 
interferogram  Figure  9* 

We  examine  the  behavior  of  the  closed  branch  of  Figure  10  made  up 
of  curves  II  and  III  for  the  case  where  |?|  4  1.  From  Eq.  (l4),  and  ^ 
the  definition  of  P,  asserting  |r|  4  1  implies  1  >  4xGWK/r  2 .  Thus  P2 
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may  be  expanded  and  only  the  first  few  terms  need  be  retained.  Eq.  (14) 
then  becomes 

-  ,  xGWK 

r  ft*  ±  7T” 
o 

The  upper  sign  refers  to  the  portion  above  the  r  =  0  axis,  the  lower 
sign  to  the  portion  below.  From  the  definitions  of  6,  and  r,  one  may 
write  for  curves  II  and  III 


2x*W2G3  K® 


(15) 


6  r  .  1  1 

|rj  W  T  W  ^  W  ’ 

for  | r|  4  1.  Hence,  Eq.  (15)  may  be  written 


(16) 


.  xGK  2x2G2K2 

-  6*77- 


(17) 


From  Eqs.  (5)  and  (8)  the  r  variation  of  a  fringe  of  the  single  fringe 
type  is  given  by 


r 


(18) 


Note  that  a  comparison  of  Eqs.  (17)  and  (l8)  shows  that  the  dependence 
on  x  of  the  first  terms  of  both  equations  is  the  same.  That  is,  for 
small  r's,  the  inner  closed  branch  should  be  similar  to  the  single 
fringe  in  appearance.  This  resemblance  is  seen  in  Figures  4  and  10. 

The  higher  order  terms  of  Eq.  (17)  cause  the  distortion  noted  in  Figure 
10,  making  the  r  values  above  the  axis  less  positive,  and  those  below 
more  negative.  There  is  some  suggestion  of  this  single-fringe -like 
closed  loop  near  the  axis  at  the  beginning  of  the  expansion  in  Figure  9. 
However,  the  processes  occuring  there  may  be  very  complex,  and  the 
simple  assumptions  resulting  in  the  closed  loop  of  Figure  10  are  very 
tenuous . 
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V.  CONCLUSIONS 


In  flow  studies  of  metal  vapor  from  a  wire  exploded  in  argon  at 
low  density,  the  effect  of  the  shocked  gas  outer  regions  can  be  neg¬ 
lected  in  the  qualitative  analysis  of  the  resulting  interferogram . 

Single  fringes  of  the  type  depicted  in  Figure  4  are  similar  to 
those  obtained  experimentally  and  shown  in  Figure  5.  Thus,  the  assump¬ 
tion  that  the  vapor  expands  in  a  hollow  cylinder  with  boundaries  that 
are  linear  functions  of  time  is  sufficient  to  explain  the  single  fringe 
behavior  observed  for  wires  exploded  in  low  density  argon.  Likewise 
the  fringes  of  Figure  7  are  very  similar  to  those  obtained  experimentally 
and  shown  in  Figure  i ,  and  it  is  sufficient  in  explaining  the  behavior 
of  the  interferograms  of  the  wires  exploded  in  a  vacuum  to  assume  a 
solid  cylindrical  expansion  with  an  outer  boundary  that  is  a  linear 
function  of  time. 

Multiple  fringes  of  the  type  seen  in  Figure  9  are  similar  to  those 
obtained  experimentally  and  shown  in  Figure  %,  lending  further  credence 
to  the  expanding  cylinder  model  for  the  explosion. 

Fringes  of  the  multiple  interferograms  may  have  closed  inner  loops 
similar  to  those  of  the  single  fringe  type,  and  Figure  9  may  show  some 
evidence  of  this  similarity. 
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Figure  1.  Flow  Assumed  in  Vaporization  Process 
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Figure  2.  Plots  of  Fringe  Shapes  Expected  from  the  Flow  Assumptions  of  Figure 
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Figure  3.  A  Plot  of  the  Function  G^) 


|G(a) 

Figure  4.  Single  Fringe  Construction. Where  a  =  0.107  and  b  =  1 
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Figure  5.  Single  Fringe  Interferogram  of  Cu  Wire 
Exploded  in  l/8  atm  Argon 


Figure  6.  Single  Fringe  Interferogram  of  Cu  Wire 

Exploded  in  a  Vacuum 


I — i 
!! 
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Figure  7.  Single  Fringe  Interferogram  for  a  =  0,  b  -  1,  and  K/6 


?-5 


Figure  8.  Multiple  Fringe  Interfere gram  Construction  for  the  r  ^  0  Case 
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Figure  9*  Multiple  Fringe  Interferogram  for  Cu  Wire  Exploded  in  l/l6  atm  Argon 


<D 
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Figure  10.  Fringe  Construction  for  rD  <  0  Cas 
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Figure  11.  Fringe  Construction  Showing  Excluded  Region 
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Figure  12.  Composite  Showing  Fringe  Constructions  Possible 

for  Multiple  Fringe  Interferograms 
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